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Abstract 

We propose an effective lagrangian for the coupling of the neutral pion with gluons 
whose strength is determined by a low energy theorem. We calculate the contribution of 
the gluonic components arising from this interaction to the pion transition form factor 
7*7* 7T° using the instanton liquid model to describe the QCD vacuum. We find that 
this contribution is large and might explain the anomalous behavior of the form factor 
at large virtuality of one of the photons, a feature which was recently discovered by the 
BaBar Collaboration. 

Pacs: 12.38.Lg, 13.40. Gp 

Keywords: pion, gluon component, instanton, form factor 



* kochelev@theor.jinr.ru 
t vicente. vento@uv.es 



The implications, of recent data by the BaBar Collaboration [I] on the transition form 
factor 7*7 — > tt°, in our understanding of the structure of the pion are being widely dis- 
cussed [21 El IH El E] • A possible scenario to explain these data consists in assuming a flat 
shape for pion distribution amplitude [31 H] , supported by some low energy models like 
the Nambu-Jona-Lasinio model [71 |8] and instanton based non-local chiral quark model 
[9], and its detailed behavior under Quantum Chromodynamics (QCD) evolution, and a 
large mass cut-off added to the quark propagator [3], signalling a peculiar behavior of the 
qq wave function [TUj. Never mind the various explanations, what appears evident is that 
the BaBar datapQ, if confirmed, are in contradiction with most model predictions based 
on the factorization approach to exclusive reactions at large momentum transfer [11] and 
this apparent lack of perturbative factorization motivates the present investigation. In 
this Letter we suggest an alternative nonperturbative explanation to the BaBar results 
based on the existence of additional contributions to the pion form factor never previ- 
ously considered. These contributions arise from the admixture of gluonic components, 
associated to nonperturbative properties of the QCD vacuum, which provide a strong 
interaction with two photons. 

Let us propose a low-energy effective ir° interaction with gluons of the following form 

C e ff = — n o^ G a Qa (1) 

J G 

Such type of Lagrangian density, describing the interaction of a pseudoscalar meson with 
gluons, was introduced many years ago by Cornwall and Soni [12J to derive Witten's 
relation between the rj' mass and the topological susceptibility, in a world without light 
quarks [13] 

f 2 
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6 

where f n = 92.3 MeV is the pion decay constant, Mx the mass of the indicated particles, 
the topological susceptibility is given by 



Xt f= =--f{M 2 r) ,+M*-2Ml\ (2) 



x ff~° = i Jd 4 x< 0|T{Q 5 (x)Q 5 (0)}|0 > G , (3) 

and 

Qb(x) = ^G%(x)G^{x) (4) 

is the topological charge density. 

The effective pion-gluon interaction, Eq{Tl is the analogue of the pion-quark effective 
interaction 

Ce J q { = ~^M q qi^fq ■ 7T, (5) 

fir 

giving the pion coupling to the quarks. 

To derive the decay constant Jq , which sets the scale of the gluon nonperturbative 
interaction with the neutral pion, we will use a low energy theorem (LET) 

'8tt ^ ^ 2m d + m t 



< 0|^G>° >= i^-^t/^. (6) 
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We stress that this matrix element is rather big due to the large light quark mass ratio 

US] 

z = ^ = 0.35 - 0.6. (7) 
By using the effective interaction EqfJJ and the LET EqE] we get 

o 2(l + z) X ?f=° 

{\-z)flMl U - W 

As it was to be expected, the strength of the coupling of the neutral pion to gluons is 
related to the violation of isospin symmetry and proportional to the difference of the d- 
and u- quark masses 

— o oc m d - m u . (9) 

JG 

For the value of the mass ratio m u jm d shown in EqJT] which is the one allowed by the 
Particle Data Group, one obtains 

R = J -f- ~ 28.1 -51.1. (10) 

J IT 




Figure 1: Gluonic contribution to the pion transition form factor. Symbol I (I) denotes 
instanton(antiinstanton) . 

Let us next calculate the contribution arising from the interaction given by EqJJJ to 
the transition form factor of two photons to the pion 7* : (qi)j* (52) _ > ^°(p), where q% and 
q2 are the photon momenta and qi + q2 = V- We consider the case where all virtualities 
of the incoming and the outcoming particles are in the Euclidean domain Q\ = —q\ > 0, 
Q\ = —q\ > 0, P 2 = —p 2 > 0. The Lagrangian density EqJT] describes the pion interaction 
with soft gluons. Such gluons should interact with photons through nonperturbative QCD 
interactions. We use the instanton liquid model (ILM) for the QCD vacuum to calculate 
the interaction of two photons with the gluonic component of the pion. The ILM is one 
of the most successful models for the description of nonperturbative QCD effects (see 
reviews [161 E])- Within the ILM the single instanton contribution to the pion form 
factor coming from the interaction Eqfljis associated to the diagrams in Fig.l. 
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The amplitude for the 7r° — 7*7* interaction via an instanton with center zq has the 
following form 



4 r 

T^{p,qi,q 2 ) = ^ap^e^—s-^ n I (p)p 2 dpx 



d 4 z J d 4 Zl J d 4 x J d^e-^e^e^QKh)) x 

^(hyA a f - h^x a ) + hr x hyx a f} , (11) 



hxhy f 1 
"A 



where 



= ^whr (12) 

is the topological charge density of the instanton in the space-time point z±, rij(p) is the 
instanton density, p is the instanton size, A = x — y, h s = l/{x 2 + p 2 ), hy = 1 /(y 2 + p 2 ) 
and the notation w = w — zq for any variable w has been introduced. The sum runs over 
the light quark flavors, i.e. i = u,d,s. To get Eq{TT]we have used the correlator of two 
electromagnetic currents in the instanton field obtained by Andrei and Gross |18j . Their 
result was corrected by a color factor (see [12]). The first term in the last line of the 
equation is coming from the quark nonzero modes in the instanton field and the last term 
arises from the interference between nonzero and zero modes [T8] . 

The final result for the gluon contribution to the pion transition form factor induced 
by instantons is 

j 4 < e 2 > 



F(P 2 , Ql QlY g = -j^- J dpnj^p'Sip, P 2 , Ql Q\) (13) 
where 

S(p, P 2 , Ql Ql) = $i(v^) f 1 dt {lit, z u Z2, z 3 ) + (1 - t)I(t, z u z 3 )} , (14) 

Jo 



Tf . \ r a a(a + l)^ 2 (Z(a,t, z u z 2 , z 3 )) 

I{t,z 1 ,z 2 ,z 3 ) = / da- — — -, 15 

Jo {a + 1 - tyZ 2 (a,t, zi, z 2 , z 3 ) 

and 

Z(a, t, zi, z 2 , z 3 ) = y{a + l)(taz 1 + tz 2 + (1 - t)z 3 )/(a + l — t). (16) 
The functions 

*i{*) = ^Y&, **{z) = zK 1 {z) (17) 

behave as $1,2 (z) — > 1 in the limit z — > 0. In Egs fTBTfTol the notations are Z\ = Q\p 2 , 
z 2 = Qlp 2 , z 3 = P 2 p 2 and < e 2 >= e 2 . 

For an estimate we use Shuryak's version of the ILM |20j, where the density is given 

by 

m(p) = n 5(p - p c ) (18) 

and 

n w 1/2/m- 4 , p c w 1/3/m. (19) 
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Within this simple model for the instanton distribution the result for the form factor is 

F{P\ Ql QlY 9 = 4< f f > rf fl S{ Pc , P\ Ql Ql), (20) 

where fj = ir 2 nop^. is so-called instanton packing fraction in the QCD vacuum. 

It should be pointed out that in spite of the smallness of instanton packing fraction 
// ~ 0.06, using the single instanton approximation as above is only valid for values 
of the momentum transfers Qi,Q2 3> 1/Rj, where Ri ~ 3p c is the distance between 
the instantons in the ILM. For smaller photon virtualities it is necessary to include the 
contributions arising from multiinstanton configurations. With an average size of the 
instanton in the QCD vacuum as in EqJT9l for the region Q\,Q\ > 1/ ' p 2 > p 2 = 0.35 
GeV 2 , i.e. zi )2 > 1, the validity of a single instanton approximation is assured. 

The calculation above was done for the case when all external momenta are Euclidean. 
In order to compare with BaBar data we have to perform an analytic continuation of the 
pion virtuality to the physical point of the pion on-shell P 2 — > —m 2 — ie. An inspection 
of the integrals in Eqs JT4]fT5l shows that the dominant contribution to the form factor at 
m 2 /Ql 2 *C 1 is coming from the region of integration t ~ 0, due to the pole at Z 2 = 0. 
Assuming the following behavior of the function § 2 (Z) = ZK\(Z) ~ 1 near the pole and 
keeping only leading terms in m 2 we obtain following closed form formulas for the real 
and imaginary parts of the flavor singlet part of form factor 



Re(F(ml,Ql,QlY g ) ~ 4< J^ /j x 



n 2 UR 

z 1 [z 2 log(z 2 )/z 1 + log(z 1 )(log(z 1 /z 2 ) - 1) + Li 2 {{zi ~ z 2 )/z 1 )} 



{ " {zi-z 2 ) 2 

z 2 \z x log{z x ) - n 2 /6 - log 2 (z 1 - z 2 )/2 - log{z 2 ){\ - log(z 1 )/2)} 



(zi ~ z 2 ) 2 

z 2 [Log(z 2 )log((zi/ (z x - z 2 )) - Li 2 (z 2 / (z 2 - Zi))} 

{zi ~ z 2 ) 2 
log(z 1 /z 2 )log(mlp 2 ^ 



z\ - z 2 



h (21) 



ImiFimlQlQDl) . ^W^ , (22) 

TXUti z x - z 2 

The imaginary part of form factor arises because the pion may decay in this calculation 
into a quark- ant iquark pair since confinement, which forbids this decay, is not explicitly 
implemented. However, the net contribution of the imaginary part to total transition 
form factor in the BaBar kinematics is very small. These formulas are useful to extract 
the behavior of the transition form factor with Q 2 . The exact numerical analysis will be 
described below. For definiteness we consider the case Z\ > z 2 

In the limit Q\ ^> Q\ which is valid for BaBar kinematics, the formulas for the real 
and the imaginary parts simplify, 

Re(F(m 7T ,Q 1 ,Q 2 )) g _ , (23) 
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It follows from Eqs j23|24l that the flavor singlet gluon induced part of the form factor 
has a dependence on the large photon virtuality Q\ proportional to log 2 (Q\) /Q\, which 
is much stronger than that of the flavor nonsinglet part, which in most of the models 
is of the form 1/Q\. The additional feature of this new contribution is its strong chiral 
enhancement since the massless logs appear governed by the pion mass as log(Q\/m 2 r ). 
For symmetric kinematics Q 2 = Q\ = Q\ the result is 



/m( ^° 2 « - Hw- (26) 

Having determined the dependence of the virtuality in the approximation Eq. |2TJ we 
proceed to study the exact numerical calculation which includes the effect of the functional 
form of ZKi(Z). Before we do so, we would like to point out that the exact calculation 
leads to a smaller (by about 50%) result, compared to the approximate calculation. This 
factor can be absorbed in the uncertainties of the vacuum model associated with the 
poor knowledge of the instanton distribution (about 30%), and the additional uncertainty 
coming from the value of the pion coupling to gluons Eqs]9] and [10] (about a factor 2) due 
indeterminacy in the ratio of u- and d- quark masses, Eqd 

We compare our result with the BaBar data. Before we do so some caveats have to 
be expressed since the comparison is not direct. The BaBar experiment, only measures 
the virtuality of one of the photons in the interval Q\ = 4 — 40 GeV 2 . They only put an 
upper limit on the virtuality for the second photon, Q\ < 0.18 GeV 2 . Finally, they use 
a model for the form factor to extract the value at the real photon point Q\ = 0. Thus, 
a direct comparison of our results with the BaBar data is not possible. Moreover, our 
calculation only represents the flavor singlet contribution to the form factor, therefore we 
have to add a flavor nonsinglet part. We take in the estimate shown in Fig. 2 for the flavor 
nonsinglet part the corresponding to a vector meson dominance (VMD) model, i.e. 

F(Ql Qir iD = * , „ 2 * , ^ 2/ ^ - (27) 



47r 2 ^(l + Q?/M 2 )(l + Q!/M 



In order to compare our results with the BaBar data we perform an extrapolation of 
their results from Q\ = to Q\ = 0.35 GeV 2 . In Fig. 2 we compare our calculation with 
the extrapolation of the BaBar data described by [TJ 

Q 2 \F°f-m\=A(-J^)\ (28) 

where A ~ 0.182 and /3 ~ 0.25. This function has been continued to the point /i 2 (0.35 
GeV 2 in our case) following the VMD model, 

f)2 I -pBaBar (r\2\ I 

Q 2 I ^-(QV) |= * \: XP 2/A , 2 )l - (29) 

The bands in the figure represent our uncertainties, both in the vacuum model and in the 
coupling constant, as mentioned before. 
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Figure 2: Contributions to the pion transition form factor compared with the extrapo- 
lation of the BaBar data to Q\ = fi 2 = 0.35 GeV 2 (solid line): gluonic contribution in 
our model including uncertainties (lower band), conventional VMD contribution (dashed 
line) line, and their sum (upper band). 



It should be mentioned that the behavior of the gluonic part of the form factor as 
function of Q 2 is determined by shape of the decay of the non-zero modes in the instanton 
field (see Eq JTT]) . At the same time it is well known that the VMD-like behavior of the 
flavor nonsinglet part of pion form factor can be easily reproduced within a non local chiral 
quark model based on the quark zero-modes dominance in the instanton vacuum pi]. Due 
to the weaker decay of the quark nonzero modes with respect to zero modes one can expect 
a harder Q 2 dependence of the flavor singlet part of the form factor in comparison with 
the flavor nonsinglet part. Such tendency is seen in Fig.2. Indeed, the gluonic part is 
described well by a fast increasing function of Q 2 , Eqj23l That function shows a similar 
behavior, as function of Q 2 , as the BaBar data, Eq]2"51 Contrary to the gluonic part of 
the form factor, its flavor nonsinglet part has a conventional 1/Q 2 behavior at large Q 2 , 
Fig.2. Taking into account some uncertainties in our estimates related to the poorly known 
ratio of the u- and d- quark masses, EqJTJ as well as uncertainties in the parameters of the 
instanton model, we conclude, that the new contribution related to the gluonic component 
of pion might explain the anomalous behavior of the pion transition form factor found by 
the BaBar Collaboration. One should be aware that these contributions are beyond the 
Operator Product Expansion and their QCD evolution is non trivial [22] . 

It is evident that such type of contribution should be present also for the i] and 
mesons. In this case one should carefully take into account the effects of their larger 
masses and their strong mixing (see recent discussion in the papers [231 121]). We also 
would like to point out that the strong Q 2 dependence of the gluonic part of the form 
factor opens a new possibility to disentangle particles with dominant gluonic content, 
i.e. glueballs (see review [25]), in 7*7* collision. These tasks are the subject of future 
investigations. 

In summary, the BaBar data [T] point towards a breaking of perturbative factorization. 
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This has led us to investigate possible non perturbative mechanism that contribute to the 
pion transition form factor. We have shown that a nonzero interaction of neutral pion 
with gluons arising from isospin violation, i.e. m u ^ m^, induces a large contribution 
to this form factor at large virtuality of one of the photons. More sophisticated models 
for the instanton density and probably multiinstanton contribution might bring the value 
closer to the observed one. 
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